home *** CD-ROM | disk | FTP | other *** search
/ The X-Philes (2nd Revision) / The X-Philes Number 1 (1995).iso / xphiles / hp48_1 / eigenvec / comp.sys.handhelds_6319_000000.msg
Internet Message Format  |  1995-03-31  |  2KB
  1. Path: en.ecn.purdue.edu!noose.ecn.purdue.edu!samsung!zaphod.mps.ohio-state.edu!sdd.hp.com!hp-pcd!hpcvra.cv.hp.com!rnews!hpcvbbs!akcs.briank
  2. From: akcs.briank@hpcvbbs.UUCP (Brian Korver)
  3. Newsgroups: comp.sys.handhelds
  4. Subject: HP48 Eigenvector program
  5. Message-ID: <281c5d39:2936.1comp.sys.handhelds;1@hpcvbbs.UUCP>
  6. Date: 29 Apr 91 18:40:06 GMT
  7. References: <sc5c=sC00WB94MH0U0@andrew.cmu.edu>
  8. Lines: 47
  9.  
  10. I think this is what you are looking for .......
  11. (User.programs) Main: read 213
  12. Item: 213 by _tasmith at hpcvbbs.UUCP
  13. Author: [Ted A Smith]
  14.   Subj: Eigenvalue/Eigenvector decomposition
  15.   Keyw: eigenvalues eigenvectors functions of a matrix
  16.   Date: Wed Feb 06 1991 22:09
  17.  Lines: 28
  18.  
  19. Here is a quick and dirty eigenvalue/eigenvector decomposition for
  20.   real symetric matricies.
  21.  
  22. I used the Jacobi method.
  23.  
  24. The termination test is a hack (I just test to see if the
  25.   eigenvector matrix has changed in a given pass!)  I don't
  26.   have any idea if there is a possibility of non-termination.
  27.  
  28. Eigen takes a real symetric matrix in level 1 and returns
  29.   the matrix of eigenvectors in level 2 and the eigenvalues
  30.   are along the diagonal of the matrix in level 1.  (The
  31.   offdiagonal values should be small in relation to the
  32.   diagonal values.)
  33.  
  34. EClr can be used to 0 the offdiagonal values.
  35.  
  36. EFun takes a real symetric matrix (M) in level 2 and a function
  37.   of 1 real arg (F) in level 1 and returns F(M) in level 1.
  38.  
  39. For example in analogy with 'SIN(x)^2+COS(x)^2==1':
  40.   [[ 1 2 3 ] [ 2 4 5 ] [ 3 5 6 ]]
  41.   DUP  \<< SIN \>> EFun DUP *
  42.   OVER \<< COS \>> EFun DUP * +
  43.  
  44.   [[ .999999999981 9.89E-12 -1.881E-11 ]
  45.    [ 1.188E-11 .999999999959 -3.3E-12 ]
  46.    [ -1.801E-11 -2.3E-12 .999999999962 ]]
  47.  
  48. ----------
  49.   Resp: 1 of 1 by _tasmith at hpcvbbs.UUCP
  50. Author: [Ted A Smith]
  51.   Date: Wed Feb 06 1991 22:11
  52.  Lines: 1
  53.  
  54. ASCII downloadable eigenvalue decomposition routines
  55.  
  56. Type attach to view and queue attached files.
  57.